A simple model for the cloud adjacency effect and the apparent bluing 

of aerosols near clouds 

Popular Summary 

In determining aerosol-cloud interactions, the properties of aerosols must be characterized 
in the vicinity of clouds. Numerous studies based on satellite observations have reported 
that aerosol optical depths increase with increasing cloud cover. Part of the increase comes 
from the humidification and consequent growth of aerosol particles in the moist cloud 
environment, but part comes from 3D cloud-radiative transfer effects on the retrieved 
aerosol properties. Often, discerning whether the observed increases in aerosol optical 
depths are artifacts or real proves difficult. The paper provides a simple model that 
quantifies the enhanced illumination of cloud-free columns in the vicinity of clouds that are 
used in the aerosol retrievals. This model is based on the assumption that the enhancement 
in the cloud-free column radiance comes from enhanced Rayleigh scattering that results 
from the presence of the nearby clouds. The enhancement in Rayleigh scattering is 
estimated using a stochastic cloud model to obtain the radiative flux reflected by broken 
clouds and comparing this flux with that obtained with the molecules in the atmosphere 
causing extinction, but no scattering. 
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29 Abstract 

30 In determining aerosol-cloud interactions, the properties of aerosols must be 

31 characterized in the vicinity of clouds. Numerous studies based on satellite observations 

32 have reported that aerosol optical depths increase with increasing cloud cover. Part of the 

33 increase comes from the humidification and consequent growth of aerosol particles in the 

34 moist cloud environment, but part comes from 3D cloud-radiative transfer effects on the 

35 retrieved aerosol properties. Often, discerning whether the observed increases in aerosol 

36 optical depths are artifacts or real proves difficult. The paper provides a simple model 

37 that quantifies the enhanced illumination of cloud-free columns in the vicinity of clouds 

38 that are used in the aerosol retrievals. This model is based on the assumption that the 

39 enhancement in the cloud-free column radiance comes from enhanced Rayleigh 

40 scattering that results from the presence of the nearby clouds. The enhancement in 

41 Rayleigh scattering is estimated using a stochastic cloud model to obtain the radiative 

42 flux reflected by broken clouds and comparing this flux with that obtained with the 

43 molecules in the atmosphere causing extinction, but no scattering. 

44 

45 



45 1. Introduction 

46 Numerous studies based on satellite observations have reported a positive 

47 correlation between cloud amount and aerosol optical thickness (AOT) (e.g., Sekiguchi et 

48 al., 2003; Loeb and Manalo-Smith, 2005; Zhang et al., 2005, Kaufman et al., 2005a, 

49 Matheson et al., 2005). Recently, Koren et al. (2007), using MODIS data, showed that 

50 the average reflectance for cloud-free ocean scenes far away from clouds were up to 30% 

51 lower than those near cloud edges. The higher reflectances lead to higher AOTs retrieved 

52 in the vicinity of clouds. This positive correlation can be explained as a result of physical 

53 phenomena such as the humidification of aerosols in the relatively moist cloud 

54 environment or a transition between aerosol and clouds where the cloud signature is weak 

55 (evaporation and/or activation of cloud drops) and the distinction between cloudy and 

56 cloud-free air becomes problematic. The term “twilight zone” was coined by Koren et al. 

57 (2007) to describe the regions around clouds which are neither precisely cloud-free nor 

58 precisely cloudy. On the other hand, part of the correlation can result from remote 

59 sensing artifacts such as cloud contamination of the cloud-free fields of view used in the 

60 aerosol retrievals. Kaufman and Koren (2006) noted that any “satellite analysis may be 

61 affected by potential cloud artifacts.” 

62 There are two ways that clouds affect the retrievals of aerosols: (i) the existence 

63 of small amounts of sub-pixel sized clouds in pixels identified as being cloud-free and (ii) 

64 an enhancement in the illumination of the cloud-free column through the reflection of 

65 sunlight by nearby clouds. When the pixels are relatively large (e.g., TOMS ~ 40 km, 

66 OMI ~ 15 km), only the first type (unresolved variability), cloud contamination is 

67 considered (e.g., Torres et al., 2002; Sinyuk et al., 2003). The second type (resolved 

68 variability), also called the ‘cloud adjacency effect,’ is more pronounced when satellite 

69 pixels are relatively small (e.g., MODIS and MISR ~ 0.5 km). Kobayashi et al. (2000), 

70 Cahalan et al. (2001), Podgomy (2003), Wen et al., (2001, 2006, 2007), Nikolaeva et al. 

71 (2005) studied the cloud adjacency effect when cloud-free pixels are brightened (or 

72 shadowed) by reflected light from surrounding clouds using 3D radiative transfer 

73 calculations applied to LANDSAT, MODIS, and ASTER data as well as to numerically 

74 generated cloud fields including an isolated cubical cloud. Both cloud contamination and 

75 the cloud adjacency effect may substantially increase reflected radiation and thus lead to 



76 significant overestimates of the AOT. These two types of cloud effects, however, have 

77 different impacts on the retrieved AOT: sub-pixel clouds increase AOT by increasing the 

78 apparent contribution due to large particles (aerosol “coarse” mode), cloud adjacency 

79 mostly increases the apparent contribution due to small particles (aerosol “fine” mode). 

80 This short paper quantifies the second factor by using a simple stochastic cloud model to 

81 obtain the radiative flux reflected by broken clouds and comparing this flux with that 

82 obtained with the molecules in the atmosphere causing extinction, but no scattering. 

83 The next section discusses the factors that contribute to the enhancement of a 

84 cloud-free column through the cloud adjacency effect. Section 3 introduces a simple 

85 two-layer model of the cloud enhancement with broken clouds as the lower layer and 

86 molecular scattering as the upper layer. A Poisson stochastic cloud model used to obtain 

87 the upward flux reflected by broken clouds is briefly described in Section 4. Section 5 

88 compares the results of this simple model with those obtained from Monte Carlo 

89 calculations for broken cumulus clouds over Brazil observed by MODIS. Finally, 

90 Section 6 summarizes the results and discusses their implications. 

91 

92 2. Cloud enhancement and its contributors 

93 Current methods used to retrieve AOT in cloud-free pixels account for sunlight 

94 reflected by the underlying surface and by the Rayleigh scattering due to molecules in the 

95 atmosphere but not the sunlight reflected by surrounding clouds. Sunlight reflected by 

96 the surrounding clouds, however, is an additional source of radiation that reaches the 

97 sensor as a result of (i) reflection by the underlying surface, (ii) scattering by the aerosol, 

98 and (iii) scattering by molecules. The relative roles of these three contributions varies 

99 from scene to scene and depends on many factors, including wavelength, surface 

100 reflectance, nearest cloud distance, cloud optical depth, the vertical and horizontal 

101 distributions of clouds, AOT, the vertical distribution of aerosols (relative to clouds), the 

102 solar and satellite viewing angles. 

103 Wen et al. (2006, 2007) gained insight into the cloud adjacency effect by 

104 performing synthesized aerosol retrievals in realistic broken cumulus fields over a 

105 biomass burning region in Brazil as observed by MODIS. They assumed that all aerosols 

106 were below the cloud tops and used 3D and ID radiative transfer calculations to 


107 determine the average difference between the 3D and ID reflectances for all cloud-free 

108 pixels as given by 

109 A p=r 3D (x,y)-r 1D . (1) 

110 The calculations were performed for a variety of surface albedos and 3 different AOTs, 

111 0.1, 0.5 and 1.0 at different wavelengths. They referred to Ap as the ‘cloud-induced 

112 enhancement’ or just ‘cloud enhancement.’ Figure 1 illustrates the results calculated for 

113 the 0.47 m wavelength. For dark surfaces the enhancement is not sensitive to AOT. 

114 For bright surfaces, the enhancement decreases with AOT because the aerosol layer 

115 prevents photons reflected by the surface from reaching the satellite. The intercept with 

116 the vertical axis gives the enhancement for zero surface albedo and thus provides 

117 estimates for the contribution from Rayleigh scattering. The contribution from molecular 

118 scattering dominates over aerosol scattering which, as is evident from the figure, is nearly 

119 an order of magnitude smaller even for an AOT of 1 . The relative roles of molecular and 

120 aerosol scattering arise because the scattering angles encountered in the retrievals of 

121 aerosol properties are typically between 100° and 150°. For this the range of angles the 

122 normalized phase functions for aerosols are much smaller than the Rayleigh phase 

123 function (e.g., Liou, 2002, p. 98). 

124 In summary, for dark surfaces and low-level clouds with aerosols below the cloud 

125 layer, sunlight reflected by the clouds and then scattered by molecules in the cloud-free 

126 columns is the key process for the enhancement of retrieved AOT, at least for the shorter 

127 wavelengths at which Rayleigh scattering is strong. Since the enhancement is due 

128 primarily to Rayleigh scattering and not very sensitive to AOT, the enhancement can be 

129 assessed knowing only the cloud properties and the (average) distance from a cloud-free 

130 pixel to a cloudy pixel. 

131 

132 3. A simple model for the cloud-induced enhancement of reflectances for nearby 

133 cloud-free columns 

134 Assume that the enhancement of the reflectance in the cloud-free column is due 

135 entirely to Rayleigh scattering. Consider a simple, two-layer model with broken clouds 

136 in the lower layer and a layer of molecules for the upper layer (Fig. 2). Take the cloud 

137 enhancement to be the difference between the following two radiances: (a) one is 



138 reflected from a broken cloud field with a scattering Rayleigh layer above it and (b) one 

139 is reflected from the same broken cloud field but with the molecules in the upper layer 

140 causing extinction, but no scattering. In other words, 


141 
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143 

v(o m r m mi a ^ d o) T Mt m (dif,e) 

l-a c (z,6 0 )R m (dif) 

(3) 

144 and 



145 

r 2 (0 o ,0) = R m (6 0 ,6) + a c (T,0 0 )T m (6 0 )t m (dif,0). 

(4) 


146 Here sub-index ‘m’ stands for ‘molecule’ while ‘c’ stands for ‘cloud.’ i? m (0o,0) is the 

147 reflectance for a molecular layer with no clouds below (this term is irrelevant here since 

148 it is canceled in calculating Ap). Cloud reflectance, 0Cc, is the critical parameter in this 

149 simple model because, in addition to cloud optical depth, T , and SZA, 0 O , it is also a 

150 function of the cloud brokenness as will be discussed below. T is the transmittance 

tn 

151 through the molecular layer with direct sunlight incident from above while t m is the 

152 transmission through the molecular layer for diffuse illumination from below. Finally, 

153 R m (dif ) is the reflectance of the molecular layer illuminated by diffuse radiation from 

154 below. Note that with the exception of a c , all the quantities in (2)-(4) are ID and are 

155 calculated using a standard plane-parallel radiative transfer code. For simplicity, the 

156 surface is assumed to be black. Contributions from non-zero surface reflectances can be 

157 readily included in a c . 

158 In summary, a simple two-layer model with a broken cloud field below and 

159 Rayleigh scattering molecular layer above is used to quantify the cloud-induced 

160 enhancement of Rayleigh scattering. The enhancement comes from the enhanced 

161 illumination of the molecular layer through the reflection of sunlight by the surrounding 

162 clouds. The main unknown is the reflectance for a broken cloud field. If we assume that 

163 the clouds are plane-parallel rather than broken then etc will be overestimated. Since Ap 

164 in (2)-(4) is an increasing function with respect to etc (Fig. 3), the plane-parallel 

165 approximation will also overestimate the effect of clouds on cloud-free pixels. 


166 The next section will describe the calculation of a c for a broken cloud field using 

167 a stochastic model. The advantage of using a stochastic model is that the output is 

168 ‘generic.’ It is averaged over many realizations of a cloud field with given statistical 

169 properties. 

170 

171 4. The Poisson stochastic model for broken clouds 

172 The one-layer Poisson model for broken clouds originally proposed by Titov 

173 (1991) is used to calculate the cloud reflectance for broken cloudy regions. Kassianov 

174 (2003) generalized this one-layer model to multilayer broken cloud fields while 

175 Zhuravleva and Marshak (2005) validated the one-layer model by comparing results with 

176 those generated using fractal cloud fields. The main parameters in the model are as 

177 follows: (0 cloud fraction, A c , (ii) averaged cloud optical depth t , and (in) cloud aspect 

178 ratio, y, which is defined as the ratio of cloud vertical to horizontal dimensions. In 

179 addition, the single scattering albedo and the cloud droplet scattering phase function 

180 along with the surface albedo are specified. For the shortwave calculations performed 

181 here, the droplet single-scattering albedo is set to unity and the Cl phase function 

182 (Deirmendjian, 1969) was used. Figure 4 shows an example of two broken cloud fields 

1 83 with A c = 0.3 and y= 0.5 and 1 . 

184 The output of the stochastic model is the domain (and ensemble) averaged upward 

185 and downward fluxes with downward fluxes subdivided into diffuse and direct 

186 components. Zhuravleva and Marshak (2005) used these subdivided fluxes to determine 

1 87 cloud aspect ratios from ground-based measurements. 

188 Note that two (averaged cloud optical depth, r , and cloud fraction, A c ) out of the 

189 three principal input parameters can be determined from the MODIS Cloud Product 

190 (MOD06). The third parameter (cloud aspect ratio y) is not readily available. 

191 Fortunately, as is shown in the next section, the cloud enhancement is not very sensitive 

192 to the aspect ratio, at least for small solar zenith angles. 

193 A simple one-layer stochastic model is used to derive cloud reflectances as a 

194 function of the average cloud optical depth, cloud fraction, and cloud aspect ratio for 

195 broken cloud regions. The clouds are distributed in space according to a Poisson 



196 distribution so that the average distance from a cloud-free pixel to a cloud edge is 

197 uniquely determined by cloud fraction and cloud aspect ratio. 

198 

199 5. Results 

200 Figure 5 shows the cloud-induced enhancement Ap as a function of cloud optical 

201 depth for 0.47 pm and four cloud fractions: A c = 1.0, 0.7, 0.5, and 0.3. The aspect ratio 

202 y = 1> the solar zenith angle 0 O = 60°, the view zenith angle 9= 0°, and the surface 

203 albedo, Os = 0.0. Note that the case of A c = 1.0 represents unbroken clouds and 

204 corresponds to the plane-parallel approximation. The figure depicts an example of a 

205 look-up-table (LUT) that can be used to estimate the expected enhancement of cloud-free 

206 radiances in the vicinity of clouds. Consider a broken cloud scene with 70% cloud cover 

207 and an average cloud optical depth of 22 illuminated by the sun with a zenith angle of 

208 60°. The enhancement in Rayleigh scattering at 0.47 pm in the nadir direction will likely 

209 be 0.04 larger than its ID counterpart. 

210 To assess the merits of the above approach, estimates of the cloud enhancement 

211 were made for the two 68 by 80 km broken cloud scenes in biomass-burning regions of 

212 Brazil studied by Wen et al. (2007). Both scenes were simultaneously observed by 

213 MODIS and ASTER. The first cloud scene (centered at 0.0N, 53.78W and acquired on 

214 Jan. 25, 2003) was described by Wen et al. (2006) while the retrieved cloud parameters 

215 for the second scene (centered at 17. IS, 42.16W and acquired on Aug. 9, 2001) were 

216 described and analyzed by Marshak et al. (2006). 

217 The first scene had cloud fraction A c = 0.53 and cloud optical depth x= 12 (std = 

218 10), and the solar zenith angle was 9 0 = 32°. The surface was covered by vegetation with 

219 a low albedo of 0.01 1 at 0.47 pm and 0.025 at 0.65 pm. For this scene, Wen et al. (2007) 

220 found an average cloud enhancement of 0.015 (std = 0.005) at 0.47 pm and 0.004 (std = 

221 0.008) at 0.65 pm (marked as ‘squares’ in the left panel of Fig. 6). Two 15 by 15 km 

222 subsets of this scene with thick (x= 14, std = 8, and A c = 0.59) and thin (x= 7, std = 6, and 

223 A c = 0.51) broken clouds were also examined using high-resolution cloud fields retrieved 

224 from ASTER data in 3D Monte Carlo simulations of the radiance fields. The cloud- 


225 induced enhancement was found to be 0.019 and 0.012 at 0.47 pm for thick and thin 

226 clouds and 0.01 and 0.0018 at 0.65 pm (marked as ‘circles’ in the left panel of Fig. 6). In 

227 addition, Fig. 6 shows asymptotic values (marked as ‘ovals’ in the left panel of Fig. 6) 

228 corresponding to the enhancements at the largest distances from cloud edges. At the 

229 greatest distances from the clouds, cloud shadows are generally avoided thereby giving a 

230 more representative estimate of the 3D effects than that obtained by averaging over all of 

231 the cloud-free pixels, some being darkened by shadows. 

232 The second scene had cloud fraction, A c = 0.4 and cloud optical depth, t = 8 (std 

233 = 8) and solar zenith angle, 9 0 = 41° (right panel of Fig. 6). The surface was much more 

234 heterogeneous than the surface for the first scene. It was also much brighter at shorter 

235 wavelengths with an average albedo of 0.04 at 0.47 pm, 0.07 at 0.65 pm (and 0.2 at 

236 0.84pm.) For this scene, Wen et al. (2007) found an asymptotic cloud enhancement of 

237 0.006 at 0.47 pm and 0.003 at 0.66pm (‘ovals’) at a distance of about 3 km from the 

238 cloud edges. The average values (‘squares’) for the cloud-free pixels selected by the 

239 MODIS AOT retrieval algorithm (Remer et al., 2005) have been included. 

240 As the results in Fig. 6 indicate, the estimates based on the stochastic model can 

241 serve as a good first-order approximation to the cloud-induced enhancement calculated 

242 with a Monte Carlo code. The stochastic model underestimates somewhat the 

243 enhancement, at least for the particular scenes studied. Clearly, the enhancement is 

244 much smaller than would be obtained with a plane-parallel approximation (A c =1). 

245 The left panel of Fig. 6 also illustrates the sensitivity of the modeled cloud 

246 enhancement, Ap, to cloud aspect ratio. For three wavelengths (0.47, 0.65, and 0.84 pm) 

247 and cloud fraction A c = 0.6 the cloud enhancement as a function of optical depth T is 

248 given for three cloud aspect ratios: y= 0.5, 1, 2. For a fixed cloud geometrical thickness 

249 of 1 km, this means that the average cloud horizontal dimension varies from 500 m to 2 

250 km. The uncertainties caused by an unknown (but reasonable) aspect ratio are of the 

251 order of 5-10%. For small cloud fractions and large solar zenith angles the modeled 

252 enhancements become more sensitive to cloud aspect ratio. 

253 The right panel of Fig. 6 also shows the effect of surface albedo. For small cloud 

254 fraction the contribution of a bright surface to the total cloud-induced enhancement can 



255 be significant. It is interesting to note that, in contrast to plane-parallel clouds, the 

256 surface contribution to the total enhancement does not decrease with cloud optical depth. 

257 It is almost constant. This is a special feature of broken cloud fields where the radiation 

258 reflected by the surface in cloud-free regions goes directly to a satellite detector rather 

259 than being attenuated by the clouds. 

260 Finally, the effect of the enhancement on the Angstrom exponent in the vicinity of 

261 clouds is studied. The Angstrom exponent characterizes the dependency of aerosol 

262 optical thickness on wavelength and is related to the average size of the particles in the 

263 aerosol: the smaller the particles, the larger the exponent. 

264 Consider three cases with the “true” Angstrom exponents equal to 0 for a “clean” 

265 environment, 1.04 for a “polluted” environment, and 2.14 for a “very polluted” 

266 environment for 0.47 and 0.65 pm, and 0, 1.31, and 2.7 for 0.65 and 0.84 pm. The clean 

267 case with zero Angstrom exponent indicates that the extinction is independent of 

268 wavelength, as it is for clouds and for large nonabsorbing aerosols, like sea salt. The 

269 AOT is taken to be 0.1 at 0.65 pm. Taking into account the cloud-induced enhancement, 

270 the “apparent” Angstrom exponent will be greater than zero. Figure 7 illustrates the 

271 increase in Angstrom exponents for both spectral intervals. Obviously, for highly 

272 polluted environments, the cloud adjacency effect is much smaller than for clean 

273 environments. Nonetheless, owing to the effects of clouds, the retrieved Angstrom 

274 exponent can be substantially larger than its true value. The cloud adjacency effect is 

275 opposite that for cloud contamination where subpixel scale clouds increase the “coarse” 

276 mode fraction thereby decreasing the Angstrom exponent. 

277 

278 6. Summary and discussion 

279 A simple model was described for estimating the cloud-induced enhanced 

280 reflectances of cloud-free columns in the vicinity of clouds. The enhancement was 

281 assumed to be due entirely to Rayleigh scattering. For the shorter wavelengths where 

282 molecular scattering is relatively large, attributing the enhancement to the illumination of 

283 the Rayleigh scattering atmosphere by sunlight reflected from nearby clouds proved 

284 reasonable (Fig. 1) for scenes with dark surfaces, broken, low-level cumulus clouds, and 


285 an aerosol layer below the cloud tops. The enhancement in Rayleigh scattering was 

286 estimated using a stochastic cloud model (Fig. 4) to obtain the radiative flux reflected by 

287 broken clouds and comparing this flux with that obtained with the molecules in the 

288 atmosphere causing extinction, but no scattering as given by (2)-(4). 

289 The results of numerical simulations of the enhancement (Wen et al. 2007) were 

290 shown to be in good agreement (Fig. 6) with the simple model, although the model 

291 underestimates somewhat the enhancement for the particular scenes studied, cumulus 

292 cloud fields retrieved from collocated MODIS and ASTER images over a biomass 

293 burning region in Brazil. 

294 The one-layer Poisson stochastic cloud model (Titov, 1990) uses cloud optical 

295 depth, r , droplet single scattering albedo and scattering phase function, cloud fraction, 

296 A c , cloud aspect ratio, y, and surface albedo to estimate reflectances for broken cloud 

297 fields. The optical depth and cloud fraction are given in the MODIS Cloud Product 

298 (MOD06). They can be used as a first approximation to quantify the cloud-induced 

299 enhancement from precalculated look-up-tables (see Fig. 5, for an example). The cloud 

300 aspect ratio is not readily available but the error due to an incorrect cloud aspect ratio is 

301 5-20%. For clouds distributed in space according to a Poisson distribution, the average 

302 distance from a cloud-free pixel to the nearest cloud is uniquely determined by cloud 

303 fraction and cloud aspect ratio. 

304 The assumption that the enhancement of the cloud-free column is due to molecular 

305 scattering leads naturally to a larger increase of AOT for shorter wavelengths, or to a 

306 “bluing” of aerosols near clouds (Fig. 7). As a result, in contrast to cloud contamination 

307 by sub-pixel clouds, the cloud adjacency effect will increase the apparent aerosol “fine” 

308 mode fraction rather than the “coarse” mode fraction. Recent findings in the MODIS 

309 cloud and aerosol products indicate that the AOT and its fine mode fraction increase in 

310 the vicinity of clouds (Kaufman et al., 2005b). 

311 Since MODIS and CERES are on the same spacecraft, another approach to 

312 estimating spectral upward fluxes for broken cloud fields is to use the CERES data. 

313 Using CERES fluxes rather than a stochastic cloud model requires the use of a theoretical 

314 radiative transfer model to convert broadband fluxes to spectral fluxes. A simpler 



315 approach would be to ignore the wavelength dependence in the anisotropy as given by the 

316 CERES Angular Distribution Models (ADMs) (Loeb et al. 2005) and use the ADMs to 

317 determine spectral fluxes from the MODIS radiances. This approach, however, can lead 

318 to the large errors at the 10 by 10 km scale of the MODIS Aerosol Product and needs 

319 further study. 

320 The enhanced illumination of cloud-free columns is a key part of characterizing 

321 aerosol properties in the vicinity of clouds. In satellite based studies of cloud-aerosol 

322 interactions, changes in the properties of the aerosol due to the cloud environment must 

323 be separated from the apparent changes that come from 3D cloud-radiative transfer 

324 effects on the retrieved aerosol properties. 

325 The simple model presented here should be taken as limited to the case of low- 

326 level clouds over dark surfaces with the aerosol below the cloud tops. The model may 

327 well prove inappropriate for scenes with highly reflecting surfaces, with upper-level 

328 clouds, or in which a substantial fraction of the aerosol lies above the low-level clouds. 

329 In such cases molecular scattering will not necessarily have the dominant role that it has 

330 for the low-level cloud and aerosol systems studied here. 

331 

332 Acknowledgments. This work was supported by the Department of Energy (under 

333 grant DE-AI02-95ER61961 to NASA’s GSFC) as part of the Atmospheric Radiation 

334 Measurement (ARM) program, and by NASA’s Radiation Program Office (under grants 

335 621-30-86 and 622-42-57). We thank C. Chiu, A. Davis, I. Koren, T. Vamai, and W. 

336 Wiscombe for stimulating discussions. 

337 


337 


References 


338 Cahalan, R. F., L. Oreopoulos, G. Wen, A. Marshak, S.-C. Tsay, and T. DeFelice, 2001: 

339 Cloud Characterization and Clear Sky Correction from Landsat 7. Remote Sens. 

340 Environ., 78 , 83-98. 

341 Deirmendjian, D., 1969: Electromagnetic Scattering on Spherical Polydispersions. 

342 Elsevier, New York (NY), 292 pp. 

343 Kassianov E., 2003: Stochastic radiative transfer in multilayer broken clouds. Part I: 

344 markovian approach. J. Quant. Spectrosc. Radiat. Transfer, 77 , 373-394. 

345 Kaufman, Y.J., L.A. Remer, D. Tanre, R.R. Li, R. Kleidman, S. Mattoo, R. Levy, T. Eck, 

346 B.N. Holben, C. Ichoku, J. Martins, and I. Koren, 2005a: A critical examination of 

347 the residual cloud contamination and diurnal sampling effects on MODIS estimates 

348 of aerosol over ocean. IEEE Trans. Geosci. Remote Sens., 43 , 2886-2897. 

349 Kaufman, Y. J., I. Koren, L. A. Remer, D. Rosenfeld, and Y. Rudich, 2005b: The effect 

350 of smoke, dust and pollution aerosol on shallow cloud development over the 

351 Atlantic ocean. Proceedings of the National Academy of Sciences, 

352 10.1 073/pnas.0505 191102. 

353 Kaufman, Y.J., and I. Koren, 2006: Smoke and pollution aerosol effect on cloud cover, 

354 Science, 313 , 655-658. 

355 Kobayashi T., K. Masuda, M. Sasaki, J. Mueller, 2000: Monte Carlo simulations of 

356 enhanced visible radiance in clear-air satellite fields of view near clouds. J. 

357 Geophys. Res., 105 , 26569-26576. 

358 Koren, I., L. A. Remer, Y. J. Kaufman, Y. Rudich, and J. V. Martins, 2007: On the 

359 twilight zone between clouds and aerosols. Geophys. Res. Lett., 34 , L08805, 

360 doi: 10. 1029/2007GL029253. 

361 Liou, K. N., 2002: An introduction to atmospheric radiation. Academic Press, New York, 

362 583 pp. 

363 Loeb, N. G., and N. Manalo-Smith, 2005: Top-of-atmosphere direct radiative effect of 

364 aerosols over global oceans from merged CERES and MODIS observations. J. 

365 Climate, 18 , 3506-3526. 

366 Loeb, N. G., S. Kato, K. Loukachine, and N. M. Smith, 2005: Angular distribution 

367 models for top-of-atmosphere radiative flux estimation from the Clouds and the 

368 Earth’s Radiant Energy System instrument on the Terra satellite. Part I: 

369 Methodology. J. Atmos. Oceanic Technol., 22, 338-351. 

370 Matheson, M.A., J.A. Coakley, Jr., and W.R. Tahnk, 2005: Aerosol and cloud property 

371 relationships for summertime stratiform clouds in the northeastern Atlantic from 

372 AVHRR observations. J. Geophys. Res., 110 , D24204, 

373 doi: 10.1 029/2005 JD006 1 65. 

374 Marshak, A., S. Platnick, T. Vamai, G. Wen, and R. F. Cahalan, 2006: Impact of 3D 

375 radiative effects on satellite retrievals of cloud droplet sizes, J. Geophys. Res., Ill, 

376 DO9207, doi:10.1029/2005JD006686. 



377 Nikolaeva, O.V., L. P. Bass, T. A. Germogenova, A. A. Kokhanovisky, V. S. Kuznetsov, 

378 and B. Mayer, 2005: The influence of neighboring clouds on the clear sky 

379 reflectance with the 3-D transport code RADUGA, J. Quant. Spectrosc. Radiat. 

380 Trans/., 94 , 405^24. 

381 Podgomy, I. A., 2003: Three-dimensional radiative interactions in a polluted broken 

382 cloud system. Geophys. Res. Lett., 30, 1771, doi:10.1029/2003GL017287. 

383 Remer, L. A., and 12 other coauthors, 2005: The MODIS aerosol algorithm, products, 

384 and validation, J. Atmos. Sci., 62, 947-973. 

385 Sekiguchi, M., T. Nakajima, K. Suzuki, K. Kawamoto, A. Higurashi, D. Rosenfeld, I. 

386 Sano, S. Mukai, 2003: A study of the direct and indirect effects of aerosols using 

387 global satellite data sets of aerosol and cloud parameters. J. Geophys. Res., 

388 108(D22), 4699, doi:10.1029/2002JD003359. 

389 Sinyuk A., O. Torres, and O. Dubovik, Imaginary refractive index of desert dust using 

390 satellite and surface observations, 2003: Geophys. Res. Letters, 30 (2), 1081, doi: 

391 10. 1029/2002GL0 16189. 

392 Titov, G. A., 1990: Statistical description of radiation transfer in clouds. J. Atmos. Sci., 

393 47 , 24-38. 

394 Torres, O., P.K. Bhartia, J.R. Herman, A. Syniuk, P. Ginoux, and B. Holben, 2002: A 

395 long term record of aerosol optical depth from TOMS observations and comparison 

396 to AERONET measurements, J. Atm. Sci., 59, 398-413. 

397 Wen, G., R. F. Cahalan, T.-S. Tsay, and L. Oreopoulos, 2001: Impact of cumulus cloud 

398 spacing on Landsat atmospheric correction and aerosol retrieval. J. Geophys. Res., 

399 106, 12129-12138. 

400 Wen, G., A. Marshak, and R. F. Cahalan, 2006: Impact of 3D clouds on clear sky 

401 reflectance and aerosol retrievals in biomass burning region of Brazil. Geosci. 

402 Remote Sens. Lett., 3, 169-172. 

403 Wen, G., A. Marshak, and R. F. Cahalan, L. A. Remer, and R. G. Kleidman, 2007: 3D 

404 aerosol-cloud radiative interaction observed in collocated MODIS and ASTER 

405 images of cumulus cloud fields. J. Geophys. Res., 23, xxxx, doi: 

406 10. 1029/2006JD008267. 

407 Zhang, J., J. S. Reid, B. N. Holben, 2005: An analysis of potential cloud artifacts in 

408 MODIS over ocean aerosol optical thickness product. Geophys. Res. Lett., 32, 

409 10. 1029/2005GL023254. 

410 Zhuravleva, T., and A. Marshak, 2005: On the validation of the Poisson model of broken 

411 clouds. Izvestiya Atmos, and Oceanic Phys ., 41 , 6, 713-725. 

412 


412 Figures 


413 

414 

415 

416 

417 



AOT=0.1 

AOT=0.5 

AOT-I.O 


Surface Albedo 


Figure 1. Cloud-induced enhancement as a function of surface albedo and AOT for a 
broken cumulus scene with cloud cover close to 50% described in Wen et al. (2007). The 
Rayleigh scattering is for 0.47 jim. 





418 Figure 2. A schematic two-layer model of a broken cloud field (lower layer) and 

419 Rayleigh scatterers (upper layer). 
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422 Figure 3. Cloud-induced enhancement as a function of cloud reflectance for four 

423 wavelengths: 0.47, 0.65, 0.84, and 1.00 p,m. The Rayleigh optical depth is taken to be 

424 0.05 at 0.65 |i.m and varies inversely with the fourth power of the wavelength. The solar 

425 zenith angle, G 0 = 60°, viewing zenith angle 0 = 0°, and the surface is black. 
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Figure 4. An example of the Poisson distribution of broken cloud fields with cloud 
fraction A c = 0.3 for a 10 by 10 km area. For a cloud vertical thickness of 1 km, the left 
panel has cloud aspect ratio y= 0.5, and the right panel has y= 1. 
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Figure 5. Cloud-induced enhancement Ap and cloud optical depth T for four cloud 
fractions, A c = 1.0, 0.7, 0.5, and 0.3. A c = 1 corresponds to the plane-parallel 
approximation. The aspect ratio is y= 1, solar zenith angle, Q q = 60°, view zenith angle, 

0 = 0°, and the surface is black. The filled circle indicates the expected cloud-free 
radiance enhancement due to nearby clouds with r -22 and A c = 0.7. 
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Figure 6. Cloud-induced enhancement Ap and cloud optical depth, t , for three 
wavelengths: 0.47, 0.65, and 0.84 pm. (Left) Cloud fraction, A c = 0.6, solar zenith angle, 
0 0 = 32°, and view zenith angle, 0 = 0°. These conditions correspond to the first broken 
Cu scene studied by Wen et al. (2007). Thick solid lines are Ap calculated using (2)-(4) 
with aspect ratio y= 1, dotted lines are with y= 2, and thin solid lines with y= 0.5. The 
surface is black. Filled blue and red squares, circles and ovals are from Wen et al. (2007) 
at 0.47 and 0.66 pm. Squares correspond to the scene average values; circles correspond 
to two subscenes with thick and thin clouds, and ovals correspond to asymptotic values. 
The dotted lines coursing through the symbols give one standard deviation. (Right) 
Cloud fraction, A c = 0.4, solar zenith angle, 0 0 = 41°, and view zenith angle, 0 = 0°. 
These conditions correspond to the second broken Cu scene studied by Wen et al. (2007) 
and by Marshak et al. (2006). The aspect ratio y= 1 . Dotted lines are Ap calculated using 
(2)-(4) for a black surface. Solid lines are for the Ap that correspond to the MODIS- 
retrieved surface spectral albedos: a s = 0.04 at 0.47 pm, a s = 0.07 at 0.65 pm, and a s = 0.2 
at 0.84pm. Filled ovals and squares are also from Wen et al. (2007). Ovals correspond 
to the actual asymptotic values while squares are the average enhancements for those 
pixels that were selected by the MODIS AOT retrieval algorithm (see text for details). 





455 Figure 7. The Angstrom exponent and cloud optical depth, t , for three situations: 

456 “clean,” “polluted,” and “very polluted.” The cloud fraction is A c = 0.5, the aspect ratio, 

457 y= 0.5, and the illumination and viewing directions are the same as in Fig. 5. Two 

458 spectral intervals are shown: left panel is for 0.47 and 0.65 |im while right panel is for 

459 0.65 and 0.84 fim. 





